201 


Appendix K \ 92-14913 

Use of Rigorous Vector Coupled-Wave Theory for Designing and 
Tolerancing Surface-Relief Diffractive Coaponanta for Magneto-Optical Haads 

Os 10 

Charles w. Haggans and Raymond K. Kostuk ✓• 

Optical Sciences Center ■') * 

University of Arizona , ' J, 

Tucson, Arizona 85721 U \ ( i 

...V 


ABSTRACT 1 

A rigorous coupled-wave model is presented, experimentally validated, and used for 
tolerancing surface-relief diffractive elements. Applications of this model in the 
design and tolerancing of components for magneto-optical (M-O) data storage heads are 
investigated. 


1. INTRODUCTION 


Surface-relief lithographically generated diffractive elements show promise .or M-0 
data Btorage head applications due to their polarization selectivity, planar geometry, 
high diffraction efficiency, and manufacturability. However, previous application of 
surface-relief diffractive elements to M-O storage systems has been limited due to the 
lack of a suitable description for their polarization properties. 


A model for describing these properties is presented in this paper. A general rigorous 
vector coupled-wave model for the calculation of the diffraction efficiency and the 
polarization properties of surface-relief gratings is described and validated with 
experimental measurements. The accuracy of the model is then analyzed. .he 
presentation concludes with examples showing how the model can be used m the design 
and tolerancing of components for M-O head applications. 


2. RIGOROUS COUPLED WAVE MODEL 


The formalism used in this model is an extension of Moharam and Gaylord's rigorous 
coupled-wave analysis for volume grating diffraction. ■*. In this treatment, an 
arbitrary surface-relief profile is approximated by slicing the profile into a number 
of equal thickness slabs. Coupled-wave equations are then generated for the field in 
each slab and are solved using a state-variables method. The complex amplitudes of the 
reflected and transmitted diffracted orders are then generated by matching boundary 
conditions . 

All fields are treated as vector quantities in this model. A vector approach enables 
the calculation of the polanation properties of diffracted beams when cross-coupling 
between the polarization states takes place. This condition exists when the gating 
vector is not in the plane of incidence. Polarization cross-coupling a 
effect that cannot be treated using scalar approaches. Vector analysis also al.ows for 
the treatment of arbitrary incident and output polarizations. 

The state of incident elliptical polarization is specified by th» s and ip amplitude 
ratio angle (a) (defined in Figure 1) and the s and p phase difference <*>•*« this 
figure, a and b are the major and minor axis lengths for the vibrational ellipse of 
incident electric vector. a, and a 2 are th* maximum field ■ t ^*ngths in P 

directions, and 0 is the angle of rotation of the major axis of the ellipse from the 
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p axis. An additional angle (*) ia used to represent the angular ratio of the manor 
to minor axes of the ellipse* These parameters are related by the trigonometric 
relations 3 : * 


cos ! A ; 


tan 2 ! V ! 

tan 2 1 <x < 


sin : A > 


ab 

a,b. 


( l>-<2) 


sin ! A i 


9in 2 ! X' 
sm 2 1 « ; 
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The output polarization parameters determined by the model are the s and p amplitude 
ratio angle (a) and the s and p phase difference angle (6) for all reflected and 
transmitted orders. From these two parameters, the rotation angle of the polarization 
ellipse (^) can be calculated. In the conical diffraction case, the s and p directions 
are different for each diffracted order. These directions are defined with respect to 
the plane containing the propagation vector for that order and the grating vector. 

The model is implemented in MATLAB, an interpreted matrix manipulation language. Runs 
are currently made on a 33MHz/386 personal computer and a VAX 8650. Typical run times 
are 3 seconds per data point for a single step grating (nt*l and ns-9) and 5 minutes 
per data point for a ten step grating (nt»10 and ns*9). 


3. EXPERIMENTAL VALIDATION OP THE MODEL 


Model predictions were compared to experimental measurements made in our laboratory and 
to the published results of two other research groups. The diffraction efficiency and 
polarization properties of a trapezoidal profile photoresist grating were measured in 
our laboratory. This grating was fabricated holographically in Shipley 1811 
photoresist. The substrate was a microscope slide. Measurements were made in 
reflection with an absorbing layer index-matched to the back surface of the substrate 
to eliminate Fresnel reflections from that surface. The incident beam was linearly 
polarized at a 45® angle to the s and p axes. The experimental reflected diffraction 
efficiency ( g ) was defined as the measured power diffracted into the *lst reflected 
order divided by the measured incident power. 

The polarization rotation angles and x were measured by both a direct and an indirect 
method. In the direct method, a linear polarizer was used to measure v. When the 
transmission axis of the polarizer was aligned with the major axis of the polarization 
ellipse, the maximum power was transmitted. For this condition, \j> was the angle 
between the polarizer transmission axis and the p axis. 

X was measured using a linear polarizer and a quarter-wave plate. The fast axis of the 
quarter-wave plate was aligned with the major axis of the ellipse. This eliminated the 
90® phase difference between the linear polarization components along the b and a axes, 
giving linearly polarized light. The plane defined by the linear polarization was 
rotated an angle x from the major axis of the ellipse. This angle was measured using 
a linear polarizer as described above. 

In the indirect method, the magnitudes of 0 and x were calculated from four power 
measurements. The measured powers were proportional to {a,;., |a.j : , !a; : , and b: : . 

! a, 1 = and |a : | 3 are proportional to the power transmitted by a linear polarizer aligned 
with the p and s axes. Similarly, powers proportional to |a| : and | b [ 2 are measured for 
a polarizer aligned with the major and minor axes. |^|, lx!, |6|, and {a; 
were then calculated using equations (l)-(5). 
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These two method* for characterizing the output elliptical polarization atate 
complement each other. The transmission axis of a single linear polarizer must be 
accurately aligned in the indirect power measurement method. In the direct 
measurements, both a linear polarizer and the fast axis of a quarter-wave plate must 
be accurately positioned. Thus, the indirect method give ,^| and ,y| more precisely 
than direct measurements. The direct measurement of and x gives the sense of 
rotation of the angles. This is useful because the sign of 6 is the same as the sign 

of x • 3 

The geometrical and physical properties of the grating were determined for input into 
the model. The grating period was determined optically to be .429 t .0005 urn. This 
value was obtained by retroref lecting the •♦•1st diffracted order onto t ^ e incident beam. 
The angle of rotation of the grating to achieve this condition (?) gives the grating 
period according to: 


In this relation, A is the grating period and X is the w * v * L * n 9 th i/llrMiaoh^sX? 
light. This value was then used to calibrate the Scanning Electron Micrograph (SEM) 
{Figure 2). The grating depth wai estimated from this photo to be .22 ± .01 um. A 1 
step approximation to the grating profile was made by averaging the profiles o< -he 
visible grating ridges and 11 coupled waves were retained in the analysis. .he index 
of * refraction of the substrate was measured to be 1.51 using a Brewster's angle 
technique, and the published Shipley photoresist index of refraction is 1.64 at .6328 

iu m. 

Figures 3 through 5 compare the experimental and calculated vaiues of diffraction 
efficiency, a, and |5| for the flat reflected order diffracted by this gracing. To 
obtain the best fit to the experimental data, iterations were mad* over the uncertainty 
rang* of the grating depth and grating profile measurements. Figure 3 *h° 

diffraction efficiency (n) versus angle of incidence. Figure 4 shows 1st order 
polarization ratio angle (a) versus angle of incidence, and Figure 5 8ho *®the phase 
difference magnitude (!«',) versus angle of incidence. Good agreement was obtained for 

all parameters. 

Our model predictions were compared to published experimental results from a study of 
the antireflection properties of short period gratings. In this work, Enger and ^ a8 
fabricated gratings etched in fused quartz substrates with periods ,hort enou f ” 
suoress all propagating diffraction orders. They observed significant phase dif .erences 
lllllln ih* s «d p polarizations of the zeroth order transmitted beam. Figure 6 show, 
the best fit to their experimental measurements and our model calculations o. phase 
dif ference J "a. t a function* of grating depth for grating, of approximately triangular 
Drofile The physical parameters of their grating no. 8 were used as input * or °^ r 

Sli Th*. S !>•. . P«tod of .31 jm. Oh. r.fr.£iv. tnd., of th. |o.r = 

substrate is 1.46, and measurements were made at a wavelength of *6328 m»® l p 

profile approximation was used and 9 coupled-waves were retained in the analysis. 

Nearly exact agreement was obtained over a 1 range of depths. 

Finally, a comparison was mad* to work originally performed by Moharam, et. al.’ and 
later by Nak?t£ and Koshiba*. Coupled wav* (C-W) and boundary-element ( B- E) theory 
were used respectively in these studies to calculate the transmitted di. fraction 
efficiency of high aspect ratio photoresist gratings. The diffraction * ° 

two gratings was measured for s and p polarized light by the original authors. The s 
and p component diffraction efficiencies for grating no. 8 wets calculated 
model. This grating has a period of .458 jim, a refractive index of 1.6 , P 

.59 jim, and an undercut profile. 

Figure 7 shows th* experimental measurements and the calculations of the three t models 
for the diffraction efficiency of . polarized light as. function 
incidsncs. Figurs 8 is ths corresponding plot for p polarisation. O 
deviated from th* measured data in th* manner fc h* previou OPO vides the 

polarization for angle, of incidence between 15- and 25«. Th* B-E model P*°^de. the 
best fit to th* measured data for large angles of incidence. However, our calculations 

0,7?.*; 3,-^L !S 

OF POO# QUALITY 



204 


Appendix K 


tand to agree with experiment better than the previous c-w model prediction*. For p 
polarization, our model provided a better fit to the experimental data for large angle* 
of incidence than either of the two previou* model*. However, it exhibited increased 
deviation* from experiment from 25® to 40®. 

Because the parameter* input to the previou* C-W model are not known, it is not 
possible to determine the cause of the poor agreement between the two model*. The 
rigorous coupled-wave formalism used in our model reduce* to the scalar coupled-wav* 
formalism used in the previous C-w model. Thus, we assume that differences in the 
predictions of the two models occur because of differences in the approximation of the 
grating profile. 


4. ACCURACY OF MODEL PREDICTIONS 


Three factors determine the accuracy of the model in predicting experimental 
measurements: The accuracy of measurement of the physical properties of the grating, 

the precision with which the grating profile is approximated, and the numerical 
precision of the model. 

For exact reproduction of the experimental measurements, the grating profile must be 
exactly characterized. Some physical parameters such as the grating period, the angle 
of incidence, and the refractive indices can be accurately determined. However, the 
grating depth and the grating profile are difficult to characterize. Measurement of 
these two parameters is limited by the irregularity of the grating ridges and by lack 
of contrast of the SEM photograph. 

The precision of the approximation of the grating profile is determined by the number 
of slabs (nt) used in the approximation. The numerical precision of the model is 
determined by the number of coupled waves (ns) retained in the analysis. For perfect 
9 r *ting characterization, each output parameter converges to a constant value as nt and 
ns are increased. 

The effects of inaccurate profile determination and insufficient numerical precision 
a *-e readily observable from Figures 9 and 10. Figure 9 shows the calculated reflected 
+ lst order diffraction efficiency versus ns for 5, 8, 10, and 12 approximations to the 
profile of Figure 2 for an incidence angle of 35®. Figure 10 is a similar plot for the 
reflected s and p phase difference. In these figures, the curves for nt*8 and nt*12 
converge to efficiency values that are larger than those for the nt*5 and nt*10 curves. 
If the shapes of the curves changed as nt increased, the precision of the approximation 
of the profile would be in question. For this situation, increasing nt would cause a 
convergence in profile shape. However, the similar shape of the curves indicates that 
error in duty cycle measurement is responsible for the different efficiencies at 
convergence. This source of error limits the accuracy of the comparison to experiment 
to approximately 5 degrees for the phase difference angle and .5% for the diffraction 
efficiency for the grating of Figure 2. 

The convergence of the curves in Figures 9 and 10 indicates that the numerical 
precision of the predictions increases as ns increases. The phase difference curves 
of Figure 10 can be observed to converge slower than the diffraction efficiency curves 
of Figure 9. This indicates that diffraction efficiency is less sensitive to the 
numerical precision of the model. Thus, the phase difference is a better indicator of 
the level of model precision. 

These figures show that the precision of the profile approximation and the numerical 
precision of the model are satisfactory for this grating for coarse profile 
approximations and a small number of retained orders. This is evident from the 
similarity of the curve shapes and the rapid convergence of each curve. Sufficient 
precision for small values of ns and nt is computationally advantageous because the 
execution time of the model increases rapidly as these parameters increase. 
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This model demonstrates that small errors in the fabrication of surface-relief gratings 
(depth, period, etc.) can cause substantial deviations from designed performance. 
These deviations are especially important as the grating period decreases and for high 
aspect ratio gratings. 

The impact of manufacturing errors can be quantified by the comparison of calculated 
diffraction efficiencies over a hypothetical range of grating parameter values. g 
11 shows theoretical first order diffraction efficiency curves for a i two step 
approximation to a trapezoidal grating for three closely spaced grating periods (1-0# 
1 P 1 and 1.2 M«t>) • This figure shows that a .2 jcm change in the grating period can cause 
U muS as a 30% change in the diffraction efficiency. Since an optimistic tolerance 
in the fabrication of grating structures using lithographic techniques is of • 
order, it can be seen that manufacturing tolerances limit the performance o. high 
spatial frequency lithographically generated gratings. 


6. APPLICATION TO POLARIZING ELEMENT DESIGN 


A uniaue aspect of this model is its three-dimensional efficiency surface generation 
in a single run. 

nonse f lective element. Figure 14 shows the indicated cross-section in the .6 grating 
period plane. 


7. CONCLUSION 


A model for calculating the diffraction efficiency and polarization properties of 
and design of diffractive components has been presented. 
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Fig. 1. Vibrational ellip»p of the electric ffolH 
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Fig. 2. ^panning electron micrograph of the 
photoresist grating measured in our laboratory. 
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Fig. 3. +lst reflected order D.E. 
for the grating in Fig. 2. 
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Fig. 7. Transmitted D.E. of grating 8 
of reference 6 for s polarization. 
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Fig.8. Transmitted D.E. of grating 8 
of refe rence 6 for p polarization. 
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Fig. 9. Reflected +1 art order DE vs. ns Fig. 10. Reflected phase difference vs. 

calculated for the grating in Figure 2. ns calculated for the grating in Fig. 2. 
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Fig. 11. +lst transmitted order diffraction efficiency for s polarization vs. angle 
of inrirtonrip for a two slab photoresist grating (n=1.64) of .55 micron depth on a 
glaaB su bs trate (n = 1.51) at a wavelength of .6328 micr ons. The duty cycles of 
the two slabs are 23.8% and 38.1% respectively, and the first slab is positioned 
on the second slab to give even symmetry in the grating profile. 




(absolute) 


Appendix K 


grating depth 


grating period 


9 - 45 * 


grating depth 


grating period 


Fig. 12. +lst transmitted order diffraction efficiency for s and 
p pola ri z ation s vs. grating depth and period for a rectangular profile 
photoresist grating (n=1.68) at a wavelength of .780 mi^mna The 
duty cycle of the grating is 60%. 
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Fig. 13. Cross-sectional profile B of 
FSg.12 for grating period=.4S microns 
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Fig. 14. Cross-sectional profile A of 
Fig. 12 for grating periods 6 microns 
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